Astronomy & Astrophysics manuscript no. Paper 


© ESO 2008 


October 5, 2008 





On the chemistry and distribution of HOC^ in M 82 
More evidence for extensive PDRs 

A. Fuente\ S. Garcia-Burillo\ A. Usero^'^, M. Gerin^ R. Neri^, A. Faure^, J. Le Bourlot^, M. Gonzalez-Garcia^, 

J.R. Rizzo^, T. Alonso-Albi^ , and J. Tennyson^ 

^ Observatorio Astronomico Nacional (OAN), Apdo. 112, E-28800 Alcala de Henares, Madrid, Spain 

^ Centre for Astrophysics Research, University of Hertfordshire, College Lane, Hatfield ALIO 9AB, United Kingdom 

^ Laboratoire d' Etude du Rayonnement et de la Matiere, UMR 8112, CNRS, Ecole Normale Superieure et Observatoire de Paris, 24 

rue Lhomond, F-75231 Paris Cedex 05, France 

^ Institut de Radioastronomie Millimetrique, 300 rue de la Piscine, F-38406 St Martin d'Heres Cedex, France 

^ Laboratoire d'Astrophysique Observatoire de Grenoble, BP 53, F-38041 Grenoble Cedex 9, France 

^ LUTH, Observatoire de Paris and Universite Paris, 7 Place Jansen, F-92190 Meudon, France 

^ Laboratorio de Astrofisica Espacial y Fisica Fundamental, Apdo 78, E-28691 Villanueva de la Canada, Madrid, Spain 

^ Department of Physics and Astronomy, University College London, Gower Street, London WCIE 6BT, United Kingdom 



ABSTRACT 

Context. The molecular gas composition in the inner 1 kpc disk of the starburst galaxy M 82 resembles that of Galactic Photon 
Dominated Regions (PDRs). In particular, large abundances of the reactive ions HOC^ and CO^ have been measured in the nucleus of 
this galaxy. Two explanations have been proposed for such large abundances : the influence of intense UV fields from massive stars, 
or a significant role of X-Rays. 

Aims. Our aim is to investigate the origin of the large abundances of reactive ions in M 82. 

Methods. We have completed our previous 30m HOC^ J=1^0 observations with the higher excitation HCO^ and HOC+ J=4^3 and 
3^2 rotational lines. In addition, we have obtained with the IRAM Plateau de Bure Interferometer (PdBI) a 4" resolution map of the 
HOC+ emission in M 82, the first ever obtained in a Galactic or extragalactic source. 

Results. Our HOC^ interferometric image shows that the emission of the HOC^ 1^0 line is mainly restricted to the nuclear disk, 
with the maxima towards the E. and W. molecular peaks. In addition, line excitation calculations imply that the HOC^ emission arises 
in dense gas (n>10'^ cm~^). Therefore, the HOC^ emission is arising in the dense PDRs embedded in the M 82 nuclear disk, rather 
than in the intercloud phase and/or wind. 

Conclusions. We have improved our previous chemical model of M 82 by (i) using the new version of the Meudon PDR code, (ii) 
updating the chemical network, and (iii) considering two diff'erent types of clouds (with diff'erent thickness) irradiated by the intense 
interstellar UV field (Go=10^ in units of the Habing field) prevailing in the nucleus of M 82. Most molecular observations (HCO^, 
HOC^, CO^, CN, HCN, HsO^) are well explained assuming that ~87% of the mass of the molecular gas is forming small clouds 
(Av=5 mag) while only ~13% of the mass is in large molecular clouds (Av=50 mag). Such small number of large molecular clouds 
suggests that M 82 is an old starburst, where star formation has almost exhausted the molecular gas reservoir. 

Key words, galaxies: individual : M 82 — galaxies: nuclei — galaxies: starburst — ISM: molecules — ISM: abundances - radio 
lines: galaxies 



1. Introduction 

M 82 is one of the nearest and brightest starburst galaxies. 
Located at a distance of 3.9 Mpc, and with a luminosity of 
3.7x10^^ Lq, it has been the subject of continuum and line obser- 
vations in virtually all wavelengths from X-rays to the radio do- 
main. Several molecular line studies indicate that the strong UV 
field has heavily influenced the physical conditions, kinematics 
and chemistry in M 82 (Mao et al. 2000, Lord et al. 1996). 

On the other hand, M 82 exhibits one of the largest opti- 
cally visible outflows or "superwinds" in the local Universe. 
Superwinds are believed to be driven by the thermal and ram 
pressure of an initially very hot (T~10^ K), high pressure 
(P/k-lO^ K cm"^) and low density wind, itself created from 
merged supernovae remnants, and to a lesser extent by the stel- 
lar winds from massive stars. This very hot component is the re- 
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sponsible for the X-ray emission detected in this starburst galaxy 
(Strickland & Heckman 2007). 

Back to the millimeter observations, M 82 is one of the most 
beautiful examples of how chemistry can help to the full under- 
standing of the interstellar medium (ISM) of an external galaxy 
and the reverse, how the extragalactic research can push forward 
the chemistry comprehension. Our early HCO interferometric 
map using the PdBI showed that the M 82 nucleus is a giant 
photon-dominated region (PDR) of -650 pc size. Furthermore 
the comparison between the HCO and H^^CO^ images sug- 
gested that the PDR chemistry is propagating across the M 82 
nucleus (Garcia-Burillo et al. 2002, hereafter Paper I). This prop- 
agation has to be understood in terms of the existence of diflTer- 
ent populations of clouds (diflferent sizes, densities and temper- 
atures) bathed by an intense UV field and with diflTerent spatial 
distribution across the M 82 disk. The densest clouds seem to be 
concentrated in two "hot spots" located ±15" from the center of 
the galaxy, E. and W. peaks (Mauersberger & Henkel 1991). Our 
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subsequent chemical study using the 30m telescope allowed us 
to put some restrictions to the ISM physical conditions. We ob- 
served the well known PDR tracers CN, HCN, HOC+, and C0+ 
towards the M 82 nucleus (Fuente et al. 2005, 2006; hereafter 
Paper II and III). The high [CN]/[HCN] ratio and the detection 
of the reactive ions HOC^ and CO^ all across the galaxy disk 
proves the exceptionally high energetic conditions prevailing in 
the ISM of this galaxy. In extragalactic research, HOC^ has only 
been detected in the active nucleus of NGC 1068 (Usero et al. 
2004) and CO^ has been only tentatively detected in Cygnus A 
(Fuente et al. 2000). 

In Paper II and III, we modeled the chemistry in the M 82 
disk using a PDR chemistry with an enhanced cosmic ray 
flux (4 10"^^ s"^). Our best model accounted for the observed 
[CN]/[HCN] and [CO+]/[HCO+] ratios in the scenario of a 
highly fragmented interstellar medium in which small dense 
cores (n~10^ cm"^, N^o^< 1.3x10^^ cm"^ ) are bathed by an in- 
tense UV field (Go =10^ Habing fields). This scenario is consis- 
tent with that proposed by Lord et al. (1996) based on obser- 
vations of FIR forbidden lines. Melo et al. (2005) catalogued a 
total of 197 young massive clusters only in the starburst core. 
This incredible density of star clusters gives rise to a very un- 
usual, highly energetic, and highly fragmented medium which is 
very diff'erent from anything in our Galaxy. The main drawback 
of our model is that it was unable to reproduce the large CO^ 
column densities measured in this galaxy. 

PDR chemistry is not the unique proposed scenario to ex- 
plain the molecular chemistry in M 82. Spaans & Meijrenik 
(2007) suggested that the high CO^ and HOC^ column densities 
could be explained if their emission arises in a X-ray dominated 
region (XDR). Recently, van der Tak et al. (2008) detected the 
molecular ion HsO^ in M 82 and concluded that the observed 
HsO^ column densities were better explained by the enhanced 
cosmic ray PDR model than by the XDR model, however. 

Our aim now is to investigate the origin, UV photon domi- 
nated or X-ray dominated chemistry, of the large abundance of 
reactive ions in M 82. In this Paper, we present a 4" angular reso- 
lution image of the HOC^ 1^0 line towards M 82, the first ever 
obtained in a Galactic or extragalactic source. In addition, we 
determine the excitation conditions of the gas emitting in HCO^ 
and HOC^. Our results discard X-rays as the major driver of the 
reactive ions chemistry in M 82. With this information in our 
hands, we improve our chemical PDR model to account for the 
molecular abundances measured in M 82 thus far. 

M 82 is one of the few galaxies in which we can resolve the 
individual molecular clouds and provides an excellent opportu- 
nity to study the PDR and/or XDR chemistries within a single 
object. This kind of studies are of paramount importance to es- 
tablish extragalactic chemical patterns to discern between the in- 
fluence of soft X-rays and UV photons on the chemistry of the 
most distant galaxies. 



2. Observations 

2.1. Single-dish observations 

Observations of the J=3^2 and J=4^3 rotational lines of 
HCO+ and HOC+ were carried out using the 30m IRAM tele- 
scope located at Pico de Veleta (Spain) and the JCMT telescope 
at the Mauna Kea (Hawaii). In Table 1 we show a summary of 
all the observations used in this Paper. The 30m telescope obser- 
vations were presented in Paper III. 
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Fig. 1. Spectra of the J=3^2 and J=4^3 rotational lines of 
HCO+ as observed with the JCMT. The W., Center and E. po- 
sitions were observed in the two lines. Intensity scale is main 
beam temperature. 

In Fig. 1 and Table 2 we show the observational results 
from the JCMT observations. All the lines were observed with a 
spectral resolution of 1 MHz. The frequencies of the rotational 
lines of both isotopomers, HCO^ and HOC^, difl'er by less than 
1 GHz. Therefore, we observed them simultaneously avoiding 
possible pointing and calibration errors in the measurement of 
the relative intensities. 

2.2. Interferomethc observations 

We present high-angular resolution observations in the contin- 
uum at 3.3 mm and in the HOC^ 1^0 rotational lines towards 
M 82. The observations were carried out with the IRAIVQ ar- 
ray at Plateau de Bure (PdBI) in 2005, May and December. The 
antennas were arranged in the C and D configurations provid- 
ing an almost circular beam of 4.4" x 3.6" (natural weight- 
ing). We adjusted the spectral correlator to give a contigu- 
ous bandwidth of 1 GHz. The frequency resolution was set to 
2.5 MHz (8.4 km/s). We subtracted the 3.3mm continuum emis- 
sion using the GILDAS software package to generate the HOC^ 
1^0 image. The noise (rms) in the spectroscopic image is of 
~1 mJy/beam. 

In addition, a 3.3 mm continuum map was generated by av- 
eraging the channels free of line emission. Uniform weighting 
was applied to the measured visibilities producing a clean beam 
of 3.5" X 2.7". The noise (rms) in the 3.3mm continuum image 
is 0.3 mJy/beam. 

The maps have not been corrected for primary beam attenu- 
ation which has at the frequency of HOC^ a HPBW of 56" . The 
images are centered at the dynamical center of the galaxy, R.A.= 
09:55:51.9, Dec=-h69:40:46.90 (J2000). 

3. Results 

3.1. Single-dish observations 

The HCO^ 3^2 and 4^3 lines have been detected towards the 
three positions observed across the M 82 disk (see Fig. 1). We 
have not detected the HOC^ 3^2 and 4^3 lines towards any 
of the observed positions using the JCMT. Gaussian fits to the 
HCO+ 3^2 and 4^3 lines and upper limits to the HOC+ lines 
are given in Table 2. 

1 IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
ION (Spain). 
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Fig. 2. a. Integrated intensity image of the HOC^ 1^0 line as observed with the PdB interferometer. Contours are 31.8, 77.7, 124.3, 
170.9, 217.6 and 264.2 mJy/beam x km/s. The inner ring is indicated by white pentagons and the outer ring by white squares, b. 
The integrated intensity map of the HOC^ 1^0 line superposed on the interferometric H^^CO^ 1^0 image reported in Paper I 
(grey scale: 100-800 mJy/beam x km/s ). c. The same as b. but on our continuum interferometric image at 3.3mm (grey scale: 1-35 
mJy/beam). d. The same as b. but on the interferometric continuum ^^CO 2^1 image reported by Weiss et al. (2001) (grey scale: 
0.7-2.5 Jy/beam x km/s). e. The same as b. but on the [Nell] image published by Achertman & Lacy (1995) (grey scale: 0.1 10-0.90 
Jy). f. The same as b. but on the HCO image published in Paper I (grey scale: 100-350 mJy/beam x km/s). 



Table 1. Line and telescope parameters 



Molecule 


Transition 


Freq. (GHz) 


Tel. 


HPBW 




J=1^0 


86.75433 


PdBI 


5.9"x5.6" 


HOC+ 


J=1^0 


89.48741 


PdBI 


4.38"x3.6" 


HOC+ 


J=1^0 


89.48741 


30m 


28" 


HCO+ 


J=3^2 


267.55753 


30m 


9" 


HOC+ 


J=3^2 


268.45109 


30m 


9" 


HCO+ 


J=3^2 


267.55753 


JCMT 


18" 


HOC+ 


J=3^2 


268.45109 


JCMT 


18" 


HCO+ 


J=4^3 


356.73413 


JCMT 


14" 


HOC+ 


J=4-^3 


357.92199 


JCMT 


14" 



To estimate the physical conditions of the emitting gas, 
mainly the molecular hydrogen density, it is essential to derive 
accurate values of the (HCO+ 4^3)/(HCO+ 3^2) line inten- 
sity (Tb) ratio (see Sect. 4 and Fig. 4). Since the two lines are 
observed with different beams, this requires of making some as- 
sumptions about the emission spatial distribution. On basis of 
our H^^CO^ 1^0 and HOC^ 1^0 interferometric images we 



have assumed that the emission of the HCO^ 3^2 and 4^3 
lines arises in a uniform slab which is unresolved in the direction 
perpendicular to the plane. In this case we only need to make a 1- 
D beam dilution correction to the intensity of the J=4^3 line to 
get the right line intensity ratio. With the assumed geometry, the 
HCO^ 3^2 line intensity ratio is -^0.8 all across the 

galaxy plane. We next discuss the quality of our approximation 
by comparing the intensities of the HCO^ 3^2 line as measured 
with the 30m telescope (HPBW~9") with those measured with 
the JCMT (HPBW- 18"). 

The linewidths of the HCO"^ 3^2 line as observed using the 
30m and the JCMT telescopes are quite similar towards the E. 
With our assumption that the emitting region is a uniform slab, 
the ratio between the line intensities should be ~2. The measured 
value is 1.4. Taking into account possible relative calibration er- 
rors (30%) between the two instruments, our measurements are 
quite compatible with the assumed geometry. 

The linewidths of the HCO^ 3^2 lines observed with the 
30m and the JCMT are, however, very different towards the 
Center. This is not surprising taking into account that the size of 
the JCMT beam is twice that of the 30m telescope (see Table 1). 



4 



A. Fuente et al.: On the chemistry and distribution of HOC^ in M 82 



Table 2. Gaussian fits to the lines observed with the JCMT 







JTmbCIv Visr 


Av 


Tm5 


Transition 


(Kkms-1) (kms-i) 


(km s-i) 


(mK) 


E. (+14",+5") 


HCO+ 


3^2 


22.6 (0.2) 296(2) 


101(3) 


210(11) 


HCO+ 


4^3 


19.9 (0.3) 305(1) 


86(2) 


217( 7) 


HOC+ 


3^2 


<1^ 


lOV 


<33^ 


HOC+ 


4^3 


<0.5 


86" 


<21 


Center (0",0") 


HCO+ 


3^2 


29.4(0.2) 200(1) 


212(2) 


130(3) 


HCO+ 


4^3 


17.6(0.3) 140(2) 


122(3) 


135(8) 


HOC+ 


3^2 


<0.4 


212^ 


<9 


HOC+ 


4^3 


<0.7 


122^ 


<24 


W. (-14",-5") 


HCO+ 


3^2 


16.8(0.3) 118(1) 


90(2) 


174(7) 


HCO+ 


4^3 


15.1(0.5) 111(1) 


73(3) 


194(17) 


HOC+ 


3^2 


<0.6 


90« 


<21 


HOC+ 


4^3 


<1.2 


73« 


<51 



^ Assumed linewidth to derive the integrated intensity upper limit. 
^ Upper limits are 3cr. 



Three diff'erent velocity components (CI centered at 255 km/s, 
C2 at 160 km/s and C3 at 95 km/s) appear in the profile of the 
HCO+ 3^2 line as observed with the JCMT (see Fig. 1). Only 
the component at 160 km/s appears in the 30m spectrum (see 
Paper III). Assuming again that the emission is uniform in a slab 
centered in the disk plane and unresolved in the perpendicular 
direction, one would expect that the intensity measured with the 
JCMT is a factor of 2 lower than that measured with the 30m be- 
cause of the diff'erent beams. The line intensity of the C2 com- 
ponent is indeed a factor -1.7 lower than the main brightness 
temperature measured with the IRAM 30m telescope. However 
this rule does not apply to the integrated intensities because of 
the contribution of the CI and C3 components in the JCMT spec- 
trum. 

Therefore, we conclude that the assumed geometry is con- 
sistent with our 30m and JCMT observations. 

3.2. Interferometric observations 

In Fig. 2 we compare our HOC^ 1^0 image with the high an- 
gular resolution images of M 82 in several molecular lines. Our 
continuum image at 3.3 mm and the [Nell] map by Achertman 
& Lacy (1995) are also shown. 

All molecular maps depict the typical nested rings morphol- 
ogy, with two maxima located almost symmetrically on both 
sides of the dynamical center of the galaxy. However, the angu- 
lar distance between both peaks is diff'erent from one species to 
another. This suggests a chemical stratification of the ring struc- 
ture. The [Ne II] emission (Achertman & Lacy 1995) and the 
radio recombination lines (Rodriguez-Rico et al. 2004) define a 
ring of HII regions with a width of ~10"-15". The peaks of this 
ring are marked with pentagons in Fig. 1 and hereafter we will 
refer to it as the inner ring. The CO and H^^CO^ maxima lies 
beyond the HII regions defining a ring with a width of ~25"- 
28". The HCO emission extends further out in the disk with a 
separation of -32" between the peaks. The peaks of this ring are 
marked with squares in Fig. 1 and hereafter we will refer to it 
as outer ring. We interpreted this chemical stratification as the 
consequence of the propagation of the PDR chemistry outwards 
the disk (Paper I). 

The emission of the HOC^ 1^0 line is weak in the inner 
ring. The HOC^ emission is concentrated in two maxima lo- 



cated at the inner edge of the outer ring, at the same distance 
from the galaxy center as the peaks of the H^^CO^ 1^0 and CO 
2^1 lines. In fact, the HOC+ 1^0, H^^CO^ 1^0 and CO 2^1 
lines have similar spatial distribution at scales of ~ 150 pc. This 
similarity breaks down at smaller scales, however. The E. max- 
imum of HOC^ is displaced ~40 pc towards the North relative 
to that of H^^CO^. Similar displacements (~40 pc) between the 
HOC"^ and H^^CO"^ maxima are observed at several positions in 
the western-half of the M 82 disk. The spatial distributions of the 
CO 2^1 and HOC^ 1^0 lines are quite diff'erent in the region 
enclosed by the inner ring. This is very likely due to the dif- 
ferent critical densities of the two transitions. Mao et al. (2000) 
claimed that the bulk of the CO emission is arising in a warm 
low density interclump medium instead of in the dense clouds. 
The spatial distribution of the HOC^ 1 ^0 line is certainly very 
diff'erent from that of the HCO F=2^ 1 fine. The HOC+ emission 
is displaced, in z-direction, towards the North and, in radial di- 
rection, towards the center of the galaxy relative to that of HCO. 
A very diff'erent spatial distribution between HCO and the other 
PDR tracers H^^CO^, HOC+, and C0+, is also observed in the 
Galactic PDR associated to the ultracompact HII region Mon R2 
(Rizzo et al. 2003, 2005). 

3.2.1. 

To quantify the observed diff'erences between the spatial dis- 
tribution of the H^^CO^ 1^0 and HOC+ 1^0 lines, we have 
calculated the map of the H^^CO^ 1^0/HOC+ 1^0 bright- 
ness temperature ratio (hereafter, R) using our PdBI images con- 
volved to the same angular resolution (5.9"x5.6"). The resulting 
image is shown in Fig. 2. R takes values between 1 and 2 in most 
of the galaxy disk without any clear trend with the distance from 
the galaxy center. This is consistent with the uniform value of R 
measured with the 30m telescope (Paper II). Larger values, R>3, 
are found in the southern half of the bubble associated with the 
supernova remnant SNR 41.95-1-57.5, which produces a North- 
South gradient in the value of R towards the western-half of the 
M 82 disk. Since the excitation conditions of the H^^CO^ 1^0 
and HOC"^ 1^0 lines are quite similar, this displacement is very 
likely due to a gradient in the [H^^CO^]/[HOC^] abundance ra- 
tio. Assuming a ^^CO/^^CO ratio of 89 and the two isotopomers 
arising in the same phase, values of R between ~1 and ~3 corre- 
sponds to [HCO+]/[HOC+] between -44 and -136. These val- 
ues are similar to those found in Galactic PDRs (Fuente et al. 
2003, Rizzo et al. 2003, Savage & Ziurys 2004) and towards the 
AGN in NGC 1068 (Usero et al. 2004) but more than one order 
of magnitude larger than those found in Galactic dense clouds 
where [HCO-']/[HOC-'] =2000-6000 (Apponi & Ziurys 1997). 



4. LVG calculations 

We have made LVG calculations to derive the physical condi- 
tions of the emitting gas. Chemical models show that high elec- 
tron abundances, X(e")>10"^, are required to have such low val- 
ues of the [HCO+]/[HOC+] ratio, <80, as measured in M 82 (see 
Paper II, Usero et al. 2004, Fuente et al. 2003). Therefore, part of 
the HCO^ and most HOC^ emission is very likely arising in the 
partially ionized cloud envelopes and/or intercloud medium, and 
collisions with electrons should be considered in our molecular 
excitation calculations. We have calculated the HCO"^ and HOC"^ 
collisional rates with electrons using the method described in 
Appendix A. While collisional rates with molecular hydrogen 
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Fig. 4. Results from our LVG calculations assuming typical physical conditions of the regions III (top), II (middle) and I (bottom) 
of the PDR (see Fig. 6 and text). In the right panels (b.,d. and f.) we show the HCO^ 4^3 /HCO^ 3^2 line ratio as a function of 
the hydrogen density (n(H2)) and the HCO+ column density, N(HCO+)= 89 x N(H^^CO+). Left panels (a., c. and e.) have to be 
interpreted as follows: the continuous line is the contour corresponding to a HCO^ 4^3 /HCO^ 3^2 line intensity ratio of 0.8. 
The HCO^ 4^3 /HCO^ 3^2 line intensity ratio does not depend on the beam filling factor because we assume that it is the same 
for the two lines. Short-dashed lines correspond to an observed H^^CO^ 1^0 intensity of 9.1 mK (value observed towards the E.) 
assuming diff'erent values of the beam filling factor for this line emission. Long-dashed lines are the contours for a HCO^ 3^2 
intensity of 210 mK (value toward the E.) assuming a [HCO^]/[H^^CO^]=89 and the same values of the beam filling factor as for 
the H^^CO^ line. Note that we only find a solution for a beam filling factor of 0.006 in region III, 0.005 assuming in region II and 
of 0.0048 in region I. These are the solutions given in Table 3. 



are of the order of 10"^^ cm"^ s"\ the collisional rates with 
electrons are ~ 10"^ cm"^ s"^ which corroborates that colli- 
sions with electrons are to be considered in our case (in which 
X(e-)>10-5). 

We have derived the physical conditions of the gas emitting 
in HCO^ by fitting the spectra derived from the H^^CO^ im- 
age smoothed to a beam of 18", and the JCMT observations. In 
our LVG analysis we have assumed a constant [HCO^]/H^^CO^] 
isotopic ratio of 89 and considered three diff'erent sets of kinetic 
temperature/electron density values that defined the typical strat- 
ification within a cloud irradiated by an intense UV field: (I) 
Ty^=1000 K, X(e-)=10-^; (II) T;t=200 K, X(e-)=10-4; and (III) 
Ty^=50 K, X(e")=0. These three regions are shown in Fig. 6 for a 
plane-parallel slab with a total visual extinction of 50 mag with 
the two sides illuminated by the intense UV field prevailing in 
M82. 

In regions II and III the hydrogen is mainly in molecular 
form and collisions with H2 dominates the excitation of the 
HCO^ molecules. However, more than 50% (the exact value de- 
pends on the exact value of the visual extinction) of the hydrogen 
is in atomic form in region I. Collisions with H could dominate 
in this region because the rates for H are usually larger than those 
for H2 (partly because H is lighter than H2). As the rates for H 
are not known for this ion, we can consider that the present para- 
H2 rates (Flower, 1999) can be used for H atoms within a factor 
of 5. Such diff'erence between H and H2 collisions rates are in- 



deed observed for CO (Flower 2001, Balakrishnan et al, 2002). 
Thus our final estimate of the molecular hydrogen density in re- 
gion I is actually a particle (H and H2) density and is accurate 
within a factor of 5. 

The results of our excitation calculations are shown in 
Table 3. Since we have observed two lines of the main isotope 
HCO^, and one of the rarer isotope, H^^CO^, we can deter- 
mine the opacity of the lines and derive the beam filling fac- 
tor, the H^^CO^ (and therefore HCO^) column density and the 
molecular hydrogen density from our LVG fitting. Our obser- 
vations are well accounted for by beam filling factors between 
0.003-0.006, which at the distance of M 82 correspond to typical 
sizes of 20-30 pc, similar to the size of a giant molecular cloud 
in our Galaxy. Note, however, that this is an area beam filling 
factor and the linear size derived in this way is only an upper 
limit to the actual cloud sizes. The H^^CO^ column density is 
N(H^^CO+)/Av~ 2 10^^ cm"^ towards the three positions. This 
value is quite robust since it is not strongly dependent on the as- 
sumed values of the kinetic temperature and electron abundance 
(see Fig. 4). The derived values of n(H2) do, however, depend 
on the assumed kinetic temperature and electron abundance. For 
negligible values of X(e") we obtain molecular hydrogen densi- 
ties larger than 5x10^ cm"^. The derived value of the molecular 
hydrogen density is lower by an order of magnitude if a kinetic 
temperature of 1000 K and an electron abundance of 10"^ are as- 
sumed (see Table 3 and Fig. 4). The derived molecular hydrogen 



6 



A. Fuente et al.: On the chemistry and distribution of HOC^ in M 82 



M 82 Cloud A^=50 mag 




- 3 



Ax (arcseconds) 

Fig.3.MapoftheHi3CO+ 1^0/HOC+ 1 ^0 brightness temper- 
ature ratio, R, (grey scale) in M 82. The galaxy has been rotated 
to make coincide its major axis with the abcissas axis. In the top 
panel we compare the H^^CO^ 1^0/HOC^ 1^0 image with 
that of the integrated intensity emission of the HOC"^ 1^0 line, 
in the bottom, with that of B^^CO^ 1^0 line (Paper I). 

Table 3. LVG results 





X(e-) 


^eff 


N(H^^CO+)/Av 


n(H2) 
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Fig. 6. Predictions for the gas temperature, electronic abun- 
dance and abundances of various species derived with our im- 
proved PDR model. The calculations have been carried out for 
a plane-parallel cloud of 50 mag, with a uniform density of 
n = 4 10^ cm"^, that is being illuminated by a UV field of 
Go = 1.0 10"^ (in units of the Habing field) from the two sides. 
The cosmic rays flux is set to ^ = 4 10"^^ s"^ We have shadowed 
the three regions discussed in Section 4 and 5. 



densities are larger than 10^ cm"^ even in the case of the extreme 
physical conditions of region I. This proves that the bulk of the 
HCO"^ emission is arising from dense gas. 

For HOC^ our information is not complete. We have de- 
tected the HOC^ 3^2 line towards the E. using the 30m tele- 
scope (Paper III). We have not detected the HOC^ 3^2 line to- 
wards any position using the JCMT. Using the 30m data and 
our interferometric HOC"^ J=1^0 image degraded to an angu- 
lar resolution of 9", we have made some LVG calculations and 
estimated the physical conditions of the gas emitting in HOC^ . 
Assuming Ty^=50 K and X(e")=0, we derive a beam averaged 
HOC+ column density , N(HOC+)/Av=2.8 lO^^ cm'^, and a 
molecular hydrogen density, n(H2)~l 10^ cm"^, towards the E. 
The derived molecular hydrogen density is only a few times 
lower than that derived from the HCO^ lines assuming the same 
physical conditions. As discussed in Section 5, these physical 
conditions are not realistic for HOC^ that is expected to mainly 
arise in region I. Assuming the most extreme physical conditions 
of region I, Ty^=1000 K and X(e")=10"^, the derived H2 density 
is ~ 2 10"^ cm"^. This large density (>10'^ cm"^) suggests that the 



HOC^ emission arises in the partially ionized cloud envelopes 
rather than in the diff'use intercloud medium proposed by Mao et 
al. (2000) to explain the CO observations. 



5. Discussion 

Following we discuss the two scenarios proposed to explain the 
large reactive ions abundances in M 82, UV-photon and X-ray 
dominated chemistry, on basis of our new observational results 
and model calculations. 



5.1. Is X-ray chemistry at work in M 82? 

Spaans and Meijerink (2007) proposed that the diff'use X-ray 
emission at 0.7 keV could contribute to produce the high abun- 
dances of reactive ions measured towards M 82. To have a deeper 
insight into the chemistry of these ions, we have compared our 
high-angular resolution HOC"^ image with the spatial distribu- 
tion of the X-ray emission. 
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Fig. 5. Comparison between the spatial distribution of the diffuse hard X-ray emission (1 .6-2.2 keV) with that of the HOC^ emission 
(left panel) and the (H^^CO^ 1^0)/(HOC^ 1^0) brightness temperature ratio (right panel). Note the HOC^ emisison is not 
spatially correlated with the diffuse hard X-ray emission nor the (H^^CO^ 1^0)/(HOC^ 1^0) brightness temperature ratio. (The 
X-ray image has been taken from Strickland & Heckman 2007). 



Strickland and Heckman (2007) reported the images of the 
continuum diffuse X-ray emission in the central 2 kpc x 2 kpc of 
M 82 for the energy bands between 0.3 to 7.0 KeV. The diffuse, 
both soft and hard, X-ray emission arises from the inner -350 pc 
(20") of the galaxy and extends along the outflow lobes. The soft 
X-ray emission (E< 1 . 1 keV) is heavily absorbed by the material 
in the galaxy plane and peaks in the outflow lobes (at z>20"). In 
contrast, the less absorbed diffuse hard X-ray emission (E< 1 . 1 
keV) is intense in the galaxy disk and extends only ~8" and ~12" 
above and below the galaxy plane. 

Since soft X-rays are heavily extinguished towards the 
galaxy plane, a detailed comparison between the spatial distri- 
bution of the soft X-ray flux and the HOC^ emission is not pos- 
sible. However, it is clear from the X-ray images that the soft 
X-rays are associated with the superwind itself and the super- 
wind sources. The superwind is generated in the inner 350 pc of 
the galaxy where intense diffuse hard X-ray emission is detected. 

In Fig. 5 we superpose the image of the diffuse hard X-rays 
emission in the band 2.2-2.8 keV published by Strickland and 
Heckman (2007) on our HOC+ and the [Hi^CO+]/[HOC+] im- 
ages. Within the disk plane, the X-ray emission is concentrated 
in the region enclosed by the inner ring where only weak HOC^ 
emission is found. In fact the two most intense HOC^ maxima 
are located just outside the X-rays emitting region. Even towards 
the center of the galaxy, the most intense HOC^ peak seems to 
avoid the peaks of the X-ray emission. 

There is no obvious relationship between the values of the 
[H^^CO+]/[HOC+] ratio and the diffuse hard X-ray emission, 
either. Unlike the X-ray emission, the [H^^CO+]/[HOC+] ratio 
does not present any systematic trend with the distance from the 
galaxy center. Moreover, the [H^^CO^]/[HOC^] ratio presents 
the highest value around SNR 41.95-h57.5, the most intense X- 
ray source of the galaxy. Therefore, morphological arguments do 



not support the scenario of X-rays being the driving agent of the 
reactive ions chemistry in M 82. 

An independent argument against the XDR scenario comes 
from the different spatial distribution of the HOC^ and SiO 
emissions in M 82. Usero et al. (2004) detected large SiO and 
HOC^ abundances towards the circumnuclear region (CND) in 
NGC 1068. In this case, the large SiO and HOC+ abundances 
were interpreted in terms of the interactions between the X- 
rays and, respectively, the dust and gas in the CND. Unlike 
NGC 1068, the SiO and the HOC+ emissions arise in differ- 
ent regions towards M 82. While the HOC^ emission arises in 
the galaxy disk (this Paper), the SiO emission is tracing the gas 
outflowing from the galaxy plane (Garcia-Burillo et al. 2001). 
This different spatial distribution suggest a different origin for 
the HOC+ and SiO emission in M 82 from that in NGC 1068. 



5.2. PDR model 

In this section, we explore the possibility of explaining the 
molecular abundances observed in M 82 in terms of PDR 
chemistry. For this aim, we have used the updated version 
of the Meudon PDR code (Le Petit et al. 2006). The chem- 
ical network has been enlarged, relative to Paper III, to in- 
clude the chemistry of HOC^. Some reactions rates have also 
been updated according to the Ohio State University database 
(http://www.physics.ohio-state.edu~eric/). In particular, the re- 
action rate of C^ -\- OH CO^ -h H has been increased to 
2.9 10"^ cm"^ s"\ which translates into a larger CO^ production 
and a better agreement between our predictions and the observa- 
tional results as we will discuss below. We have also included 
a new C0+ formation reaction, O -h CH+ C0+ -h H, but its 
contribution to the total CO^ formation is very low (<3%). 



8 



A. Fuente et al.: On the chemistry and distribution of HOC^ in M 82 



In Paper III, we fitted the [CO+]/[HCO+], [HCO+]/[HOC-'] 
and [CN]/[HCN] ratios with a plane-parallel model, assuming 
00=10"^ in units of the Habing field, n=4 10^ cm~^, a cosmic ray 
flux of 4 10"^^ s"^ and clouds with a total thickness of ~ 13 mag. 
The thickness of the clouds was estimated as twice the thick- 
ness of the plane-parallel one side illuminated PDR that fits the 
observations. Our model was quite succesful in predicting the 
[CO+]/[HCO+], [HCO+]/[HOC+] and [CN]/[HCN] abundance 
ratios but could not account for the measured CO^ column den- 
sities. We needed about -20-40 clouds along the line of sight to 
explain the CN column densities, but then the predicted CO^ col- 
umn densities fell short by one order of magnitude. Conversely, 
if we fitted the CO"^ column densities, the predicted CN column 
densities were far in excess. This option looked less plausible 
because of the large number of clouds required to fit the CO^ 
column densities (>100). 

With our new chemical network, we obtain CO^ column 
densities a factor of ~5 larger than before (see Table 5). To im- 
prove our model, we have assumed , in addition, a plane-parallel 
cloud illuminated from the two sides. This is a more realistic 
view of the ISM in M 82. The parameters of our PDR model are 
shown in Table 4. 

In Fig. 6 we show our model predictions for a plane-parallel 
cloud with a total visual extinction of 50 mag. The physical 
and chemical conditions are symmetric since the cloud is illu- 
minated from the two sides by the same UV field. We can distin- 
guish three regions across the cloud according to their chem- 
ical properties: (I) Ay =0.0 1-0. 5 mag from the cloud surface. 
CO"", HCO+, HOC"", CN and HCN present an abundance peak in 
this region and [CN]/[HCN]~20 . (II) Av=0.5-5 mag from the 
cloud surface. The abundances of C0+, HCO+, HOC+ decline 
in this region. CN and HCN are, however, very abundant and 
[CN]/[HCN]~100. (Ill) Av>5 mag from the cloud surface. The 
reactive ions CO"^ and HOC"^ have negligible abundances but 
HCO+, CN and HCN are abundant with [CN]/[HCN]~1. Note 
that the reactive ions CO^ and HOC^ only survive in region I. 
Therefore, reactive ions are excellent tracers of the number of 
clouds since their total column density remains constant regard- 
less of the thickness of the cloud for clouds with Av>2 mag. 

In Tables 5, 6 and 7, we compare the predicted molecular 
column densities and abundance ratios with those observed to- 
wards M 82. In addition to C0+, HOC+, HCO+, CN, HCN and 
HCO, we include the HsO^ column density derived by van der 
Tak et al. (2008) and the upper limit to the H2O column den- 
sity obtained by Wilson et al. (2007). Thus, the most relevant 
molecules for the PDR/XDR chemistry are considered. 

In Table 6 we show a comparison between our one-type- 
cloud model and the column densities measured towards the E. 
position in M 82 (Paper I, II, III; van der Tak et al. 2008). This 
is the position towards which the emission of most molecules 
peak and column densities are better determined. The number 
of clouds along the line of sight have been calculated using the 
HOC^ column densities. To match the observed HOC^ column 
density, we use 44 clouds of Av=5. This choice provides a good 
fit (within a factor of 3) to many of the other observed column 
densities (see Table 7). However, the model falls short by 2 or- 
ders of magnitude to predict the HCN column densities. Related 
to it, the [CN]/[HCN] ratio is -380, far in excess the observed 
value. The model does not reproduce the HCO and H2O column 
densities either. Moving to ~10 mag clouds, we need 40 clouds 
along the line of sight to fit the HOC^ column density. In this 
case, we improve our fit to the HCN column density and the 
[CN]/[HCN] ratio, but the predicted CN and H3O+ column den- 
sities are then too large. This is the same problem we found in 



Paper III. We cannot fit the CN, HCN, and reactive ions column 
densities simultaneously. The new reaction rates improve the fit 
but not enough. The fit would become worse if we go to larger 
clouds. The [CN]/[HCN] ratio would decrease and get closer to 
the observed value but the predicted CN, HCN, HCO+ and H3O+ 
column densities would increase and become much larger than 
the observed values. 

In order to look for a better solution, we have tried a two- 
cloud-type model. We obtain a good agreement with the observa- 
tions by assuming that we have ~44 clouds, 98.5% of which have 
a thickness ~5 mag but a few percentage, 1.5%, have a thick- 
ness of 50 mag. Since the assumed cloud morphology (a col- 
umn plane-parallel cloud) is not realistic, the number of clouds 
has here only a statistically meaning. It means that the emis- 
sion can be explained with a cloud distribution such that most of 
the molecular gas (87%)is forming small clouds (represented by 
our plane-paralell Av=5 mag cloud) while only a few percentage 
(13%) is in larger molecular clouds (represented by our plane- 
paralell Av=50 mag cloud). Compared with the case of 100% 
of small clouds, the small percentage of large clouds inject the 
lacking HCN column density and push down the CN/HCN ra- 
tio to fit the observed value of ~6 without changing that much 
the [CO+]/[HCO+] and [HCO+]/[HOC+] ratios. In addition, our 
predicted HCO column densities are much closer, although still 
below, the observed value. Note that the number of large clouds 
is limited by the HCN and HCO^ column densities. Both species 
are very abundant in region III of our pattern 50 mag cloud and 
their column densities increase rapidly with the number of large 
clouds. With the proposed mix, we fit all the observed column 
densities except those of HCO and H2O. 

Our model falls short by a factor of 10 to predict the ob- 
served HCO column densities. Schilke et al. (1992) proposed 
that the photodissociation of the H2CO molecules released from 
the grain mantles could increase the HCO abundance in the mod- 
erately illuminated regions of the PDR. Our model does not in- 
clude production of these species on grains. 

Regarding H2O, the predicted column densitites are 2 orders 
of magnitude larger than the upper limit obtained by Wilson et 
al. (2007). This discrepancy is very likely because we are us- 
ing a gas-phase model that neglects the depletion of water in 
dense clouds. It is also important to recall that this limit was 
obtained with a large beam (2.V) and by observing the 556.936 
GHz transition of ortho-H20 which is very likely optically thick. 
Therefore, the limit is quite uncertain. 

Summarizing, we conclude that the most relevant PDR 
molecules are well fitted by assuming that most of the mass of 
the molecular gas (87%) is forming small clouds with a total 
thickness of ~5 mag. A small percentage, however, is encom- 
passed in large molecular clouds, -50 mag, in which star for- 
mation might be proceeding. This suggests that the mass distri- 
bution of the clouds in M 82 is very diff'erent from that in our 
Galaxy, and most of the giant molecular clouds (GMCs) have al- 
ready been destroyed by the recent starburst. M 82 is, therefore, 
an old starburst where star formation has almost exhausted the 
molecular gas reservoir. 



6. Conclusions 

In this paper, we present new JCMT observations of the HCO^ 
and HOC"^ 3^2 and 4^3 lines and the PdB interferometric im- 
age of the HOC+ 1^0 line towards the nucleus of M 82. Our 
results can be summarized as follows: 
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Table 4. PDR model 
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Table 5. Model Results (I) 
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- We have carried out excitation calculations for HCO^ and 
HOC"^ taking into account our 30m and JCMT observations 
of the HCO+ and HOC+ J=3^2 and J=4^3 lines. Our re- 
sults show that the two ions arise in dense gas (n>10'^ cm"^). 
Therefore, HOC^ is likely tracing the partly ionized and 
moderately dense envelopes of molecular clumps, rather than 
in the intercloud medium or wind phase. 

- We present a high angular resolution (4") image of HOC^ 
1 ^0 in M 82, the first ever obtained in Galactic and extra- 
galactic research. The comparison between the diff'use X-ray 
emission and our HOC^ image discards X-rays as the ma- 
jor driver of the high abundances of reactive ions in M 82. 
UV photons from massive stars seem, therefore, the major 
drivers of the molecular chemistry in M 82. 

- Our PDR chemical model shows that most molecular abun- 
dances measured towards M 82 (HCO+, HOC+, C0+, CN, 
HCN, HsO^) are well accounted in the scenario of two-types 
(diff'erent thickness. Ay =5 mag and Ay =50 mag) of dense 
clouds (> 10"^ cm"^) bathed by an intense interstellar UV field 
(Go=104). 

- Furthermore, our model shows that most of the molecular 
gas, 87 % in mass, is forming small clouds (Av~5 mag) while 
only 13% is in large molecular clouds (Av~50 mag). This 
suggests that only a small percentage (<13%) of the molec- 
ular gas in M 82 is regions where massive star formation is 
still proceeding. 
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Appendix A: Rotational rate coefficients 

For the rotational excitation of HCO"^ by hydrogen molecules, 
we employed the rate coefficients computed by Flower (1999). 
These rates are available for transitions between levels with 
/ < 20 and temperatures in the range 10 < T < 400 K. For 
HOC^, de-excitation rates were assumed identical to those of 
HCO"^ and the detailed balance principle was applied to derive 
the HOC^ excitation rates. Note that these rates correspond to H2 
in its ground para level (/ = 0) and that rates for H2 in excited 
rotational levels are not available. 

Rate coefficients for electron-impact rotational excitation of 
HCO"^ and HOC"^ were computed using the /^-matrix method 
combined with the Adiabatic-Nuclei-rotation (ANR) approxi- 
mation. The ab initio /^-matrix model of Faure & Tennyson 
(2001) for HCO"^ was employed here for all calculations on both 
molecular ions. The (theoretical) equilibrium geometry of HOC^ 
was taken from Ma et al. (1992). A ground state dipole moment 
of 2.797 D was obtained, in excellent agreement with the value 
of 2.8 D computed by Defrees et al. (1982). The two lowest tar- 
get states of symmetry and were included in the close- 
coupling calculation. Fixed-nuclei T-matrices were obtained for 
total symmetry ^S^, and ^A, as in the case of HCO^. Full 
details about the 7?-matrix model will be published elsewhere. 

The major difference with the calculations of Faure & 
Tennyson (2001) concerns the application of the ANR approxi- 
mation. This latter has been recently shown to be valid down to 
threshold for molecular ions, provided a simple "Heaviside cor- 
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Table 6. Model Results (II) 
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Table 7. Model Results (III) 

M82 100% A, =5 mag 100% A,=10mag 98.5% A,=5 mag 
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* Average value in our HCO interferometric image (Garcia-Burillo et al. 2002) 

** We mark in bold face those predictions that differ by more than a factor of 3 from the observed values 

rection" is applied to excitation cross sections. Full details can 
be found in Faure et al. (2006). Note that closed-channel effects 
are neglected here as these are expected to be small for strongly 
polar ions (see Faure et al. 2006). Furthermore, in contrast to the 
calculations of Faure & Tennyson (2001) restricted to the lowest 
3 rotational levels, we have considered here rotational transitions 
between levels with / < 20 (i.e. involving rotational energies up 
to ~900 K). Transitions were however restricted to A J < 8 owing 
to the limited number of partial waves included in the T-matrices 
(/ < 4). Finally, de-excitation rates have been fitted over the ana- 
lytic form, Eq. (2), of Faure et al. (2004) for temperatures in the 
range 5 < T < 1000 K. Fitting coefficients for both HCO+ and 
HOC"^ will be published in a forthcoming publication (Faure & 
Tennyson, in preparation). Note, finally, that the present HCO^ 
rates have also been employed recently by Jimenez-Serra et al. 
(2006) to explain the overexcitation of HCO^ in the shock pre- 
cursor component of the young LI 448 -mm outflow. 



